Dissipative preparation of entangled states between two spatially separated 

nitrogen-vacancy centers 
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We present a novel scheme for the generation of entangled states of two spatially separated 
nitrogen-vacancy (NV) centers with two whispering-gallery-mode (WGM) microresonators, which 
are coupled either by an optical fiber-taper waveguide, or by the evanescent fields of the WGM. We 
show that, the steady state of the two NV centers can be steered into a singlet-like state through a 
dissipative quantum dynamical process, where the cavity decay plays a positive role and can help 
drive the system to the target state. The protocol may open up promising perspectives for quantum 
communication and computation with this solid-state cavity quantum electrodynamic system. 

PACS numbers: 03.67.Bg, 42.50.Pq, 78.67.-n 



I. INTRODUCTION 

Quantum entanglement is among the most fascinat- 
ing aspects of quantum mechanics, and entangled states 
of matter are now widely used for fundamental tests of 
quantum theory and applications in quantum informa- 
tion science [ij. Plenty of different systems have been 
investigated to faithfully and controUably prepare entan- 
gled states of matter, among which cavity QED d i] 
offers one of the most promising and qualified candidates 
for quantum state engineering and quantum information 
processing. The general approach for quantum state en- 
gineering in cavity QED is based on either unitary dy- 
namical evolution [4|-ll3| or dissipative quantum dynam- 
ical process |14l - [26| . For the latter approach, the inter- 
action of the system with the environment is employed 
such that dissipation drives the system into the desired 
state. This process can be achieved by engineering the 
dissipative dynamics such that the desired state is the 
only steady state regardless of the initial state. 

Recently, the solid-state counterpart of cavity QED 
system consisting of NV centers in diamond and WGM 
microresonator has attracted great interests l28| . 
which circumvents the complexity of trapping single 
atoms and can potentially enable scalable device fabri- 
cations. This composite system takes the advantage of 
both sides of NV centers and WGM microresonators, i.e., 
the exceptional spin properties of NV centers [1^ and 
the ultrahigh quality factor and small mode volume of 
WGM microresonators [sol - lssj . Therefore, the applica- 
tion of this solid state cavity QED system in quantum 
state engineering and quan tum information processing 
is of great interests [3J-[41[. In particular, it has inter- 
esting applications in quantum networking and quantum 
communication since NV centers coupled to WGM 
microresonators are the natural candidates for quantum 
nodes, and these nodes can be connected by quantum 
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channels such as optical fiber-taper waveguide. 

In this work, we propose an efficient scheme for the 
preparation of entangled states of two spatially separated 
NV centers with two microsphere resonators (MRs) cou- 
pled either by an optical fiber-taper waveguide, or via 
the evanescent fields of the WGM. This proposal ac- 
tively exploits the resonator decay to drive the system 
to a singlet-like entangled stationary state through an 
engineered reservoir for the NV centers. We show that, 
the steady state of the two NV centers can be steered 
into a singlet-like state through a dissipative quantum 
dynamical process. Up to our knowledge, this is the first 
proposal for preparing entanglement of distant NV cen- 
ters employing reservoir engineering. Compared to pre- 
vious works using unitary dynamical evolution (36l . [37| , 
the present work has the following distinct features: (i) 
it performs well starting from an arbitrary initial state, 
which renders the initialization of the system in a pure 
state unnecessary; (ii) the produced entangled state is 
a pure stationary state; this auspicious feature is very 
promising in view of the quest for viable, long-lived en- 
tanglement; (iii) since cavity decay has been used to ac- 
tively drive the system dynamics, it thus converts a detri- 
mental source of noise into a resource. This work may 
represent promising steps toward the realization of en- 
tanglement with the solid-state cavity QED system. 



II. DISSIPATIVE ENTANGLED STATE 
PREPARATION VIA RESERVOIR 
ENGINEERING 

A. Two WGM microresonators connected by an 
optical fiber-taper waveguide 

We first consider two negatively charged NV centers 
coupled to two remote MRs connected by a fiber-taper 
waveguide, as shown in Fig. 1. NV centers in diamond 
consist of a substitutional nitrogen atom and an adjacent 
vacancy having trapped an additional electron, whose 
electronic ground state has a spin S' = 1 and is labeled 



2 



as ^^2) = l^^o) Irris = 0, ±1), where \Eo) is the orbital 
state with zero angular momentum projection along the 
NV axis. Quantum information is encoded in the spin 
states \ms — ±1) of the '^^2 triplet such that |0) = \ms — 
— 1), and |1) = |tos = +1). The A three-level system 
could be realized in the NV center if the excited state |e) 
is chosen as |e) = -^{\E-)\ms = +1) + |£'+)|ms — -1)) 

[IE 113 1 where \E±) are orbital states with angular mo- 
mentum projection ±1 along the NV axis. 
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FIG. 1. (Color online) (a) The schematic of two NV centers 
coupled to two MRs, respectively, which are connected by a 
fiber-taper waveguide, (b) Energy level structure with cou- 
plings to the cavity mode and driving laser fields. Quantum 
information is encoded in the spin states \ms — ±1) of the 



•^A2 triplet, i.e., |0) 
excited state is |e) - 



V2 



rus = — 1), 
i\E-)\ms 



and |1) = \ms = +1). The 
= +1) + \E+)\ms = -1)). 



The modes of spherical resonators can be classified by 
mode numbers n, / and m, which determine the charac- 
teristic radial (n) and angular {I and m) field distribution 
of the modes. Usually the so-called fundamental WGM 
(n = 1,1 = m) attracts great interests, whose field is con- 
centrated in the vicinity of the equatorial plane of the 
sphere. In perfect spheres, resonance frequencies do not 
depend on m, which means the MR supports both clock- 
wise and counter-clockwise fundamental WGM. However, 
it has been shown that this degeneracy can be lifted ei- 
ther by shape effect or by internal backscattering [4^ , 
and these modes can be selectively excited by coupling 
to a tapered fiber (soj . Therefore, in this work we only 
consider a single fundamental WGM with frequency 1^0 
interacting with the NV centers. 

The fundamental WGM dispersively couples the tran- 
sition |1) o |e) for each N-V center with coupling con- 



stant gj. The laser field of frequency wq couples the 
transition |0) O |e) with Rabi frequency flj. These two 
transitions are assumed to be on the two-photon Raman 
resonance, thus establishing Raman transitions between 
the states |0) and |1) through the WGM and laser field. 
The laser fields of frequencies uji and lu2 drive the tran- 
sitions |0) O |e) and |1) O |e) with Rabi frequencies 
Aj and Ilj, which are also on Raman resonance. There- 
fore, another Raman laser system is established between 
the states |0) and |1) involving the two laser fields with 
frequencies Wi and uj2. The transition |1) o |e) is also 
driven by a third classical laser field of frequency with 
Rabi frequency Ej. This laser field is used to induce an 
extra Stark shift for the state which can break the 
symmetry of the system Hamiltonian to ensure that the 
system has a unique steady state 18]. The corresponding 
detunings for the related transitions are Aj — Ueo — uq — 

UJel - Vo, ^'j = l^eO -i^l= Wei " ^2, (5j = Wei " W3, whcrC 

Weo, Wei are the transition frequencies for the NV centers. 
In the interaction picture, the Hamiltonian describing the 
interaction between the NV centers and the laser fields 
and cavity modes is (let h— 1) 

^1= ^[l],|e),(0|e»^^*-|-A,|e),(0|e^^> 

+H,|e),(l|e^^;*-fI],|e),(l|e*^'^*+5,a,|e),(l|e^^^*]-HH.c., 

(1) 

where aj is the annihilation operator for the ]th cavity 
mode. 

We now discuss the coupling between the two MRs 
and the optical fiber-taper waveguide. We consider an 
optical fiber-taper waveguide near the equatorial planes 
of both microspheres [3l| . with length I and the decay 
rate of the resonator' fields into the continuum of the 
fiber modes v. The number of longitudinal modes of the 
fiber that significantly interact with the corres pon ding 
resonator modes is on the order of Iv/2itc (qI. [lQi.Tr3. Il3| . 
If we consider the short fiber limit Iv jl-Kc < 1, then only 
one resonant mode b of the fiber interacts with the cavity 
modes. Then the interaction Hamiltonian describing the 
coupling between the cavity modes and the fiber mode is 



J^^j = i/6(a| -I- e'^fal) + H.c, 



(2) 



where v is the cavity-fiber coupling strength, and ip is 
the phase due to propagation of the field through the 
fiber. Define three normal bosonic modes c, ci,C2 by 
the canonical transformations c — "^(o-i ~ e^*'''a2), ci = 

^{di+e'''fa2 + V2b),C2 = i(ai + e-*'^a2-\/2S). In terms 
of the bosonic modes c, Ci and 62, the interaction Hamil- 
tonian J^cj is diagonal. We rewrite this Hamiltonian as 
J^o = y^iyciCi — \/2i'C2C2. So the whole Hamiltonian in 
the interaction picture is ^ = Jffo + J^i . 

In the following, the system-environment interaction is 
assumed Markovian, and then is described by a master 
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equation in Lindblad form 
dp 



dt 



with 



i[p, M'\ + &R,p+ ^R^p + %p + ^sponP, (3) 

SrjP = Kj{2djpa!j - a!jdjp - pd]aj), 
S^fp^ Kf{2bpP ~b^bp~ pPh), (4) 

where Kj is the leakage rate of photons from resonator 
j {ki — K2 — k is assumed in the foUowing) and k/ 
is the decay rate of the fiber mode. The term ^sponP 
describes spontaneous emission of the NV centers from 
the excited state |e). Its concrete form is not relevant 
to this proposal, since we will adiabatically eliminate the 
excited state in the following. 

We proceed to perform the unitary transformation 
e"'*^"*, which leads to 



+ A,|e),(0|e^^* 



+ i5i|e)i(l|e'^i*(cie-'^''* + cae*^'^* + V2c) 

+ ^52|e)2(l|e^^^*(cie-'^''* + cae*^"* - V2c) + H.c. 



(5) 



If we assume v ^ {Aj, A^, Sj^gj}, we can safely neglect 
the nonresonant modes ci , C2 under the rotating- wave ap- 
proximation. Thus we obtain the following Hamiltonian 



iAit 



A,| 



3 = 1,2 

+n,|e),(l|e*'^^*-fE,|e),(l|e 
1 



.Sit] 



(6) 

In the large detuning hmit, {|Aj|, |A^|, |(5j|, |Aj — 
A;.|, |A, - <5,|, \S, - Ar\} » {|^|, |A,|, |n,-|, \g,\}, we can 
adiabatically eliminate the excited state |e), and get an 
effective Hamiltonian describing two distinct Raman ex- 
citations 



canonical transformations b = -^{ci — C2),ai = i(ci -|- 
C2 + \/2c), a2 = ^e'^'^{ci + £2 — \/2c), we have 

!^Rj^p = K{2dipa[ — d\dip — pa[di) 

= |(ci -I- C2 + V2c)p{c\ + 4 + V2c') 

-^(ci + 4 + V2c^)(ci + C2 + V2c)p 

-^p(c| + 4 + V2c^)(ci+C2 + \/2c). (8) 

Since under the rotating-wave approximation the non- 
resonant modes ci , £2 are mostly in the vacuum state 
and have been neglected, we can discard the terms 
Cipc!j,clcjp, pc\cj,Cipc\cpc\,c\cp, c^Cip, pc\c, pc^Ci. The 
same considerations can be applied to f^/fjP and S>fp. In 
this case, S^Rip + &r2P + ^fp can be approximated as 

M 



^rP — K{2cpc^ 



c^cp - 



pc^c) 



Now, the master equation Q reduces to 
dp 



-^=i[p,Jf] + &RP 



(9) 



(10) 



We can introduce the photon number representation for 
the density operator p with respect to the normal mode 
c, i.e., p = J2'^^ri=oPmn\m){n\, where pmn are the field- 
matrix elements in the basis of the photon number states 
of the normal mode c, and are still operators with re- 
spect to the NV centers. We now assume that the res- 
onator modes are strongly damped, in which case the 
populations of the highly excited modes can be neglected. 
Therefore, we consider only the matrix elements pmn in- 
side the subspace {|0), |1)} of the photon numbers. Un- 
der this approximation, the master equation (|10p leads 
to the following set of coupled equations of motion for 
the field-matrix elements [i^ 

Poo = Fpoa - «.9cffPoi|l)i(0| - «goffPoi|l)2(0| 

+j5off|0)i(l|/Oio + i.9cff|0)2(l|pio + K/5ii,(lla) 

pii = Fpii - igeffPio|0)i(l| - j5cffPio|0)2(l| 
+igcff|l)i(0|poi + «3cff|l)2(0|poi - Kpii, 

(lib) 



J^ = -[Ai|l)i(l|+A2|l)2(l|]-[e|l)i(0| 

+ e|l)2(0|+5eff'C^|l)l(0|-hgcffC^|l)2( 



H.c. 



(7) 



whereA, = ^,e= AlHl^AlHl 



^3 - S, ' ^ - A[ - ' ana ffoff - y2Ai ~ 

^/^A ' derivation of Hamiltonian ([7]) , the resonance 



condition 



lA.I^ 



2A,- 



ic'c) is used, and 



A,- ' A' — A' 

j J J ^ 

the constant terms has been omitted. 

We now consider the dissipation term &R^p + &R^p + 
^fP, which can be simplified if we give out its expression 
in terms of the bosonic modes c, ci and £2 . From the 



Poi = Fpm - igoffPoo|0)i(l| - j.gcff/Ooo|0)2(l| 

+«5cff|0) 1(1 1/911 + «gcff|0)2(l|/Oii - -^poi, 



(11c) 

where 

Fpij = -i[M'd,Ptj\ 

= ie[|i)i(o| + |o)i(i| + 11)2(01 + 10)2(11, p.,] 

+i[Ai|l)i(l|+A2|l)2(l|,P.,] (12) 

The reduced density operator for the NV centers is 
qn = Tri?(p) ~ Poo+Pii- For the case of strong resonator 
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damping, we can adiabatically eliminate the elements poi 
and pii from the above equations [4^, prompting the 
evolution of the two NV centers with an effective master 
equation 



dg 



N 



dt 



i[gN, y^d] + ^eQN^e 



eQN 



where 



2^ 



-[Ai|l>i(l| 

+e|o)i(i|4 



+ A2|i)2(i| + e|i)i( 
e|o)2(i|] 



(|l>i(0| + |l>2(0|). 



(13) 

+ 0|1)2(O| 

(14) 
(15) 



The first term in the right side of master equation (|13p de- 
scribes the coherent laser driving of the two NV centers, 
while the last two terms describe an effective engineered 
reservoir for the NV centers. 

In order to gain more insight into the combined ef- 
fect of the unitary and dissipative dynamics described 
by the master equation ([T3| . we switch to the collec- 
tive state picture for the two NV centers, with the three 
triplet states {|00) = |0)i|0)2, |11) = |l)i|l)2,|r) = 
1/a/2(|0)i|1)2 + |l)i|0)2)}, and the singlet state \S) = 
1/a/2(|0)i|1)2 - |l)i|0)2). In this case, the Hamiltonian 
M'd and the Lindblad operator can be rewritten as 



^d - -[%/2e|ii)(r| + V2e|r)(oo| - a|^)(t| + H.(dp) 



?cff 



(|ii)(r| + |r)(oo|). 



(17) 



where we have taken Ai = — A2 = A. The effective pro- 
cesses described by the Hamiltonian ,^d and the Lind- 
blad operator is shown in Fig. 2. The Hamiltonian 
corresponding to coherent laser driving induces tran- 
sitions between the three triplet states {|00), |11), 
and between the triplet state |T) and the singlet state 
IS"), while the Lindblad operator will drive the tran- 
sition from |00) to |T), and from |T) to |11). Hence, the 
combined effect of the unitary and dissipative dynamics 
drives essentially much of the population to the singlet 
state 1 5), and a minor overlap with |11). It can readily 
be seen that the state 



\^s) = 



1 



v/2e2 + A2 



= (A|11) +\/2e|^)) (18) 



is the unique stationary state 
(HI i.e., 



ISj of the master equation 



Qjqit ^ 00) = 1-05) (V'sl 



(19) 



Therefore, starting from an arbitrary initial state, we can 
prepare the steady state of the two NV centers in the 
singlet-like entangled state |?As)- H should be noted that 
a similar master equation to (|13p has been presented in 
Refs. [H, [1^, where they study the preparation of a 
maximally entangled state of two trapped atoms in a 




FIG. 2. (Color online) EflFective processes of the system. 
The driving causes transitions between the three triplet 
states. The NV centers decay through the cavities from |00) 
to |T), and from |r) to |11) with an effective decay rate 
Kcff ~ 9,g1g/K.. The singlet IS) is coherently coupled to \T) 
by the level shift A. 



heavily damped cavity. As pointed in Refs. [3 [13 ; 
to guarantee that the state jV's) is the unique stationary 
state of the master equation (1131) , breaking the symmetry 
in the system accomplished by a small energy level shift 
is necessary. 



Two coupled WGM microresonators through 
the evanescent fields of the WGM 



We point out that the essential idea of above discus- 
sions can be applied to the system in which the two WGM 
resonators couple with each other directly via coherent 
photon hopping [i^ll^], i.e., coupled resonators. It can 
also be applied to the situation where two negatively 
charged N-V centers are positioned near the surface of 
a common WGM resonator ^6] . In the following we will 
discuss the former case in more detail. 



Photon hopping 



NV1 O 




0NV2 
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FIG. 3. (Color online) The schematic of two NV centers 
coupled to two MRs, respectively, which are coupled via the 
evanescent fields of the WGM 

We now investigate the setup in which the two WGM 
resonators interact with each other directly via the 
evanescent fields of the WGM rather than an optical 
fiber-taper waveguide, as illustrated in Fig. 3. The cou- 
plings between the resonator modes and the laser fields 
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are just the same as that shown in Fig. 1(b), except 
that the detuning for the WGM is now taken as A^. In 
this case, the interaction between the NV centers and 
the laser fields and cavity modes is just described by the 
Hamiltonian ([T]). The Hamiltonian describing the direct 
coupling between these two resonators is 

J^Q^ J{a\a2 + aial), (20) 

where J is the hopping rate of photons between the res- 
onators. At present, the master equation describing the 
system-environment interaction is 

^ = ^] + &Rip + &R2P + ^sponP, (21) 

where = + ^1. 

Introducing two normal modes di = l/\/2(ai + ^2), 
and 6,2 = l/-\/2(ai — 02), then the Hamiltonian will 
have the diagonal form — J{d\di—d}2d2)- Performing 
the unitary transformation e~*'^o* will lead to 

+n,|e),(l|e^'^^* + I],|e),(l|e^^^*] 

+ -^gi|e)i(l|e'^i*(die-*-^*-Kd2e^-'*) 

+ ^g2|e)2(l|e'^^*(rfie-*''* - ^e^-^*) + H.c(22) 
v2 

Taking A,- J = A„ |A,+ J| » {|A,|, |A;.|, |<5,|,5,}, this 

condition enables us to neglect the nonresonant mode d2 
under the rotating-wave approximation. Consequently 
we obtain the following Hamiltonian 

-fH,|e),(l|e'^^* + S,|e),(l|e^^^*] 

+ -i=gi|e)i(l|e^^i*di + -Lg2|e)2(l|e'^^*di + H.c. 

(23) 

If we further assume {|Aj |, |A' |, |Aj - A^.|,|Aj - 
5jl \5j - A;-|} > |Aj|, \Yiji then we can adi- 

abatically eliminate the excited state |e), and get an ef- 
fective Hamiltonian in a rotating frame as that of Eq. 


= -[Ai|i)i(i| + A2|i)2(i|] - [e|i)i(o| 

+e|l)2(0| +5eff4|l)l(0| +.9cffdl|l)2(0| +H.C.] 

(24) 

Following the same reasoning as that for getting the mas- 
ter equation (jl3p , we can get an effective master equation 
which has the same form as Eq. (IT5|) . In this case, we can 
also prepare the stationary state of the two NV centers 
in an entangled state. 



III. NUMERICAL SIMULATIONS AND 
PRACTICAL CONSIDERATIONS 

It is necessary to verify the model through numeri- 
cal simulations. To provide an example, we consider the 
tapered-fiber case. The same results can be obtained in 
the case where the microresonators are directly coupled. 
Taking account of the dephasing effect, we simulate the 
dynamics of the system by the following master equation 

^ = i[p, + SiR^p + SIr^p + &fp + ^dephP, (25) 
where 

2 

^dophP = 7v ^{'^^lP(^i - ^WzP - P^Wz) (26) 

with the pure dephasing rate of the NV centers, and 
ai = |l)j(l| — |0)j(0|. In the following through solving 
the full master equation (|25[) numerically, we calculate 
the time evolution of the fidelity, and of the population 
of the singlet state and of the three triplet states for 
different initial state. The fidelity with respect to the 
singlet-like state \^s) is defined as F = (V'slewlV's)- In 
order to perform the simulation much easier, we switch to 
the field representation with respect to the normal modes 
c, ci , C2 . Then the Hamiltonian in equation (|25p is taken 
as ([5]) . To solve the master equation numerically, we em- 
ploy the Monte Carlo wave function formalism from the 
quantum trajectory method ^48'|. All the simulations are 
performed under 50 trajectories. The relevant parame- 
ters are chosen such that they are within the parameter 
range for which the scheme is valid. 




time (1/g) 



FIG. 4. (Color online) Time evolution of the fidelity F, of 
the respective population of the singlet state and the three 
triplet states Ps, Pars, Pii, Pt from two different initial states, 
(a) jOO); (b) |11). The parameters are chosen as g\ ^ g2 — 

g,u ~ lOOg, Aj A^ ~ Wg,5i ~ -S2 ~ 20g,ni ~ -^2 ~ 

g, Aj ~ g, Ilj ~ O.lg, ~ VEg/5, k = 0.5g, 7^ = 0. 

Figure 4 shows the numerical results for the fidelity F, 
the population of the singlet state Ps , and the respective 
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population of the three triplet states Pqq, Pu, Pt for an 
appropriate set of parameters, starting from two different 
initial states |00) and without dephasing. One can 
readily see that the model performs very well with the 
chosen parameters. Starting from a given initial state, 
the populations Pqo , Pt oscillate rapidly with an envelop 
decaying at a rate proportional to g1^j/K, while the fi- 
delity F and the populations Ps and Ai converge to the 
maximum values at the same rate. In fact, the system 
can evolve to the stationary state jV's) with a fidelity 
higher than 99% regardless of the initial state. Figure 5 
displays time evolution of the fidelity under several de- 
phasing rates starting from two different initial states 
1 00) and |11). From the simulations we can see that, 
when 7^ < {Q^gl^j/n}, the effect of dephasing on the 
fidelity of the scheme can be neglected. However, as the 
dephasing rate 7^ becomes comparable to the effective 
coherent coupling strength 0, or larger than it, the fi- 
delity of the scheme is spoiled by dephasing. Therefore, 
to implement this proposal with high fidelity requires 
that 7,^ < {Q,glff/n}. 




o.a\ ^ , 1 

600 1200 

time (1/g) 

FIG. 5. (Color online) Time evolution of the fidelity F un- 
der several dephasing rates starting from two different initial 
states, (a) jOO); (b) |11). Other parameters are the same as 
those in Fig. 4 

For experimental implementation of this protocol, the 
present day achievements in the experiment with NV cen- 
ters and WGM resonators can be used. Coherent cou- 
pling between individual NV center in diamond and the 
WGM in microsphere or microdisk resonator has been 
reached [13, HI]. For example, in the experiment demon- 
strating the normal mode splitting with NV centers in 



diamond nanocrystals and silica microspheres [27'], the 
coupling strength between NV centers and the WGM 
is about ff/27r ~ 55 MHz, and the cavity decay rate is 
k/2ti ^ 50 MHz. For the NV centers, the electron spin 
relaxation time Ti of diamond NV centers ranges from 
several milliseconds at room temperature to seconds at 
cryogenic temperature. For our scheme, it should be im- 
plemented at cryogenic temperatures ranging from 6 to 
12 K [l^]. The dephasing time T2 induced by the fluctu- 
ations in the nuclear spin bath has the value of several 
microseconds in general, which can be increased to 2 mil- 
liseconds in ultrapure diamond [2^. Therefore, the time 
for reaching the stationary state ]V's)j which is about 
1000/(7 ^ 3/^8, is shorter than the typical decoherence 
time for this system. The two resonators can be coupled 
either by an optical fiber-taper waveguide for the case 
of remote resonators, or via the evanescent fields of the 
WGM for the case of nearby resonators. As for the coher- 
ent coupling between the fiber and resonators, a perfect 
fiber-cavity coupling (with efficiency larger than 99.9%) 
can be realized by fiber-taper coupling to high-Q silica 
microspheres [sol ] or microtoroidal cavities |49l ]. Thus 
with present day techniques in solid-state cavity QED, 
this proposal could be realized. 

IV. CONCLUSION 

We have presented a novel scheme for the preparation 
of entangled states of two spatially separated NV centers 
with two WGM resonators coupled either by an optical 
fiber-taper waveguide, or via the evanescent fields of the 
WGM. The proposal actively exploits the resonator decay 
to drive the system to a singlet-like entangled stationary 
state through an engineered reservoir for the NV centers. 
This work may represent promising steps toward the real- 
ization of entanglement with the solid-state cavity QED 
system. 
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